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The developmental role of the T-box transcription factor Tbx1 is exquisitely dosage-sensitive. In this study, 
we performed a microarray-based transcriptome analysis of E9.5 embryo tissues across a previously gener- 
ated Tbxl mouse allelic series. This analysis identified several genes whose expression was affected by Tbxl 
dosage. Interestingly, we found that the expression of the gene encoding the cardiogenic transcription factor 
Mef2c was negatively correlated to Tbxl dosage. In vivo data revealed Mef2c up-regulation in the second 
heart field (SHF) of Tbx1 null mutant embryos compared with wild-type littermates at E9.5. Conversely, 
Mef2c expression was decreased in the SHF and in somites of Tbx1 gain-of-function mutants. These results 
are consistent with the described role of Tbx1 in suppressing cardiac progenitor cell differentiation and in- 
dicate also a negative effect of Tbx1 on Mef2c during skeletal muscle differentiation. We show that Tbx1 oc- 
cupies conserved regulatory regions of the Mef2c locus, suggesting a direct effect on Mef2c transcription. 
However, we also show that Tbx1 interferes with the Gata4^ Mef2c regulatory pathway. Overall, our study 
uncovered a target of Tbx1 with critical developmental roles, so highlighting the power of the dosage gradi- 
ent approach that we used. 



INTRODUCTION 

The T-box transcription factor TBX1 is encoded by the main 
gene haploinsufficient in DiGeorge syndrome, which is char- 
acterized by congenital heart defects, hypo/aplasia of the para- 
thyroid and thymus glands, craniofacial dysmorphism as well 
as learning and behavioral abnormalities (1-5). Mutation of 
murine Tbxl can model the DiGeorge syndrome phenotype 
(6-8). A progressive dosage reduction in Tbxl mRNA is asso- 
ciated with the non-linear increase in phenotypic severity (9) 
while over-expression of Tbxl results in structural heart and 
thymic defects (10), confirming that Tbxl function, during 
embryonic development, is exquisitely dosage-sensitive. 
Uncovering the genetic and phenotypic changes resulting 
from Tbxl haploinsufficiency might help dissection of 
molecular mechanisms underlying the DiGeorge syndrome 
etiology. 



Loss of Tbxl is associated with reduced proliferation and 
premature differentiation in the second heart field (SHF) 
(11-14). Accordingly, FgfS and FgflO were down-regulated 
in Tbxl~'~ mutants (13,15,16) while genes required for myo- 
cardial differentiation (like Raldh2, Gata4 and TbxS) were 
ectopically expressed (13). 

The mechanisms by which Tbxl regulates its targets are not 
all clear; there are now examples of transcription-dependent 
and transcription-independent actions. For example, in vitro 
and in vivo evidences show that Tbxl positively regulates 
expression of the Vegfr3 gene through interaction with a 
T-box-binding element (TBE) (17), as expected for T-box 
transcription factors. 

On the other hand, it is also clear that Tbxl can function in a 
transcription-independent manner. Indeed, Tbxl is able to 
negatively regulate the bone morphogenic protein-Smadl 
pathway by binding Smadl and thus preventing the Smadl- 
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Smad4 interaction (18), which is required for the activation of 
Smadl target genes (19). 

Furthermore, Tbxl can interact with the serum response 
factor (Srf) and promote its proteasome -mediated degradation 
(17), which, in turn, results in cardiac actin and a-smooth 
muscle actin protein down-regulation (17,20,21). This non- 
transcriptional interaction, although not yet completely clari- 
fied, probably contributes to the inhibition of cardiomyocyte 
differentiation observed in vivo (13,14). 

In this study, starting from microarray-based gene expres- 
sion data, collected across an allelic series of Tbxl embryo 
mutants, we found that Tbxl negatively modulates Mef2c 
through several mechanisms. 

Mef2c belongs to the MEF2 (myocyte-specific enhancer- 
binding factor 2) subfamily of MADS [MCM1 (Minichromo- 
some Maintenance 1 Protein), AG (Agamous), DEFA (Defi- 
ciens) and SRF (Serum Response Factor)] transcription 
factors (22). Mef2 genes are expressed in cardiac, smooth 
and skeletal muscle cells, endothelial cells as well as in a 
restricted set of other tissues (23-25). Targeted inactivation 
of the Mef2c gene in the mouse resulted in cardiac and vascu- 
lar defects and embryonic lethality at E9.5 (26-28). In par- 
ticular, in Mef2c~ f ~ mice, the heart tube does not undergo 
looping morphogenesis, the future right ventricle does not 
form and several cardiac muscle genes are down-regulated. 
These data indicate a key role of Mef2c in the transcriptional 
pathways controlling myocyte differentiation. 

We provide in vivo and in vitro data, indicating that Tbxl 
plays an inhibitory role onto the Gata4^Mef2c regulatory 
pathway, and we propose that this is a mechanism by which 
Tbxl affects muscle cell differentiation. 

Our data also highlight the power of multiallelic gene ex- 
pression analysis as a tool to identify developmentally critical 
genes that correlate genotype-dependent phenotypic changes 
to genome-wide transcriptional features. 



RESULTS 

Identification of Tbxl dosage-sensitive transcripts 

We have previously generated a series of genotypes associated 
with a nearly continuous variation of Tbxl mRNA dosage 
between 0 and 1 00% of the wild-type (wt) level by combining 
two different hypomorphic alleles, Tbxr eo2 and Tbxl Neo , and 
a null allele, Tbxl~ (9). To identify Tbxl dosage-dependent 
genes in vivo, we performed a microarray-based transcriptome 
analysis across the allelic series; RNA isolated from whole 
E9.5 embryos with nine different genotypes (9) was hybri- 
dized to Affymetrix GeneChip Mouse Genome 430 2.0 
arrays. For each genotype, we used two embryos, each hybri- 
dized to one array; thus, we analyzed a total of 1 8 arrays (the 
whole set of results has been deposited in NCBI's Gene Ex- 
pression Omnibus and is accessible through GEO Series ac- 
cession number GSE33064; http://www.ncbi.nlm.nih.gov/ 
geo/ query/ acc.cgi?token=rxiphqgaeyioipy&acc=GSE3 3 064) . 

We developed a data analysis strategy to identify genes that 
respond significantly to dosage variation of Tbxl mRNA (see 
Materials and Methods) and we found 2230 transcripts with 
significant differential expression (Supplementary Material, 
Table SI). At a second filtering step (see Materials and 



Methods), we narrowed down the list to 497 transcripts (in 
addition to Tbxl) differentially expressed across the allelic 
series (Supplementary Material, Table S2). Among these 
genes, we found very limited overlaps with previously 
reported gene lists derived from Tbxl mutants (13,29); these 
genes are highlighted in green color in Supplementary Mater- 
ial, Table S2. However, those studies were carried out using 
different developmental stages and different experimental 
approaches that did not use allelic series. 

We applied a Gene Ontology (GO) search to the more strin- 
gent data set using DAVID (30,31) and found 14 GO categor- 
ies that are significantly enriched (Supplementary Material, 
Table S3). 

We also carried out a cluster analysis that returned four 
main clusters (Fig. 1). Cluster 2 includes Tbxl and transcripts 
that have an expression trend similar to Tbxl across geno- 
types. Cluster 3 behaves in an opposite way as transcripts in 
this group tend to be more expressed as Tbxl expression 
goes down. Clusters 1 and 4, which are also the richest in tran- 
scripts (see Supplementary Material, Table S4, for transcript 
content of clusters), are very different and have somewhat 
complementary patterns. In these two clusters, there appears 
to be a threshold-like effect at dosages between experimental 
points 4 and 5, corresponding to genotypes Tbxl +, ~ and 
jt) X l Neo2,Ne ° 2 > which are associated with a small overall differ- 
ence in Tbxl RNA dosage (10-15%) (9), but genetically, they 
have an important difference, i.e. the former has an intact wt 
allele and the latter has none. How this genetic difference 
can translate into the expression patterns represented in Clus- 
ters 1 and 4 remains to be investigated. 

Tbxl negatively regulates Mef2c expression 

Among the 497 differentially expressed genes, we noted that 
Mef2c expression increased as Tbxl dosage decreased (simi- 
larly to Cluster 3 of Fig. 1; Supplementary Material, 
Table S2). Mef2c is important for heart development and car- 
diomyocyte differentiation. The heart of Mef2c homozygous 
mutant mice does not undergo looping and has gross abnor- 
malities in the outflow tract, while the right ventricle fails to 
form (28), suggesting a key role for MeJ2c in SHF develop- 
ment. In particular, it regulates transcriptional pathways con- 
trolling differentiation. So, the negative modulation of Mef2c 
expression by Tbxl could contribute to explain the role of 
Tbxl in suppressing cell differentiation (14). 

We confirmed Mef2c up-regulation in Tbxl ~'~ embryos by 
quantitative real-time reverse transcriptase -polymerase chain 
reaction (qRT-PCR) of mRNA from E9.5 whole embryos 
(Fig. 2A). In addition, in situ hybridization (ISH) on embryo 
cryosections at E9.5 revealed up-regulation and ectopic 
expression of Mef2c in the SHF of Tbxl~'~ mutants 
(Fig. 2F-I) compared with wt littermates (Fig. 2B-E). If 
Mef2c expression is Tbxl dosage-dependent in vivo, we 
expected to observe a Mef2c down-regulation in Tbxl 
gain-of-function mutants. Whole-mount ISH was performed 
on CQfcT;Mespl Cre/+ embryos at E9.5; the COET (for Condi- 
tional Over-Expression of Tbxl) transgene expresses Tbxl 
mRNA in response to Cre recombination (10). Mespl Cre 
drives Tbxl expression in mesodermal domains of the 
pharynx, head mesenchyme and splanchnic mesoderm 
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Figure 1. Cluster analysis of differentially expressed genes across the allelic series. Lines represent the average profile of each cluster and the shaded areas 
represent the error (SD). Profiles were obtained by plotting the log ratios of all dosage points with respect to the control (wt) point for all genes, and then, 
the average profile, for each cluster, was determined. Tbxl profile in Cluster 2 is highlighted by black circles. The transcript content of each cluster is listed 
on Supplementary Material, Table S4. 



(Fig. 3A-D). Whole-mount ISH showed a strong reduction in 
Mef2c somite expression in CQET;Mespl Cre/+ embryos 
(arrowheads in Fig. 3E and H) and a mild Mef2c down- 
regulation in the SHF (ovals in Fig. 3G and L) or splanchnic 
mesoderm (arrows in Fig. 3E, H, F and I). This different 
effect upon Mef2c expression in different tissues may be due 
to biological reasons (different mechanisms of regulation in dif- 
ferent tissues) or to technical reasons (e.g. different levels of ex- 
pression of the COET transgene or different efficiency of 
Cre-mediated excision). We confirmed that Mef2c expression 
is Tbxl dosage-dependent using qRT-PCR of mRNA from 
C2C12 myoblasts transfected with different amounts of a 
7&x7-expressing vector (pCDNA-Tbxl-c-myc) (Fig. 4A and 
B). Western blotting (WB) analysis of nuclear protein extracts 
from transfected C2C12 cells with an anti-Mef2c antibody con- 
firmed negative regulation at a protein level, too (Fig. 4C). 



expression of myogenic markers MyoD and myogenin at day 
2 of differentiation (data not shown). Also Mef2c expression 
was increased at day 2 and, interestingly, transfected Tbxl 
inhibited this up-regulation (Fig. 4D). These data are consist- 
ent with those obtained in vivo with the gain-of-function 
model (Fig. 3E-L) and suggest that Tbxl may inhibit skeletal 
muscle differentiation through MeJ2c down-regulation. To 
further support this hypothesis, we transfected C2C12 cells 
with Tbxl, let them differentiate for 5 days and test different 
skeletal muscle differentiation markers by qRT-PCR. 
Results showed that while Myogenin did not change in re- 
sponse to Tbxl transfection, all the genes encoding muscle 
structural proteins tested (Myhl, 2, 8 and Actal) were signifi- 
cantly down-regulated (Supplementary Material, Fig. SI). 
However, myotube formation did not appear to be grossly 
affected (data not shown). 



Tbxl interferes with Mef2c expression and muscle cell 
differentiation 

We tested the response of Mef2c expression to Tbxl dosage 
during skeletal muscle differentiation. C2C12 myoblasts can 
be induced to differentiate, in vitro, by switching culture con- 
ditions from a medium supplemented with 10% fetal bovine 
serum to a medium supplemented with 2% horse serum 
(32). Under these conditions, we observed an increased 



Tbxl inhibits Gata4- but not Isll-induced up-regulation of 
Mef2c activity 

To test the effects of Tbxl upon Mef2c transcriptional activity, 
we performed luciferase assays using a MEF2 luciferase report- 
er construct (fifes , MEF2-luc) on transfected C2C12 cells (33,34). 
With this reporter, luciferase activity is proportional to active 
MEF2 transcription factors. Mef2b is expressed at a very low 
level in these cells (data not shown); Mef2a and Mef2d are 
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Figure 2. Meftc is up-regulated in Tbxl loss-of-function mutants. (A) qRT-PCR analysis showed Meftc up-regulation in Tbxl-mx\\ mutant embryos compared 
with wt littermates at E9.5. Values are means + standard error (SE), n = 3, *P < 0.05. (B-I) Meflc is up-regulated (arrows) and ectopically expressed (arrow- 
heads) in the SHF of Tbxl mutants (G and I) compared with wt littermates at E9.5 (C and E), as shown by section ISH. (C and E) and (G and I) are high- 
magnification views of the same wt (B and D) and mutant sections (F and H) respectively. 



Tbxl Mef2c 




Figure 3. Meflc is down-regulated in Tbxl gain-of-function mutants. (A-L) Left- and right-side views of whole wt (A, B and E-G) and COET; Mespl " + 
mutants (C, D, H, I and L) at E9.5. Tbxl is over-expressed in the mesodermal domains of the pharynx, head mesenchyme and splanchnic mesoderm of COET; 
Mespl Cre/+ mutants (C and D) compared with wt littermates (A and B). Meflc is mildly down-regulated in the SHF (ovals in G and L) and splanchnic mesoderm 
(arrows in E, F, F£ and I), and more extensively in the somites (arrowheads in E and H) of COET; Mespl Cre + mutants (H, I and L) compared with wt littermates 
(E, F and G). 
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Figure 4. Tbxl negatively regulates Mef2c expression in C2C12 myoblasts. (A-C) Mef2c mRNA and Mef2c protein levels are affected by Tbxl dosage in 
transfected C2C12 myoblasts as confirmed by qRT-PCR (A and B) and WB analysis (C). (D) qRT-PCR analysis showed Tbxl -dependent inhibition of 
Mef2c up-regulation during in vitro C2C12 myoblast differentiation. Values are means + SE, n = 3, *P < 0.05. 



detectable, but their expression is not significantly affected by 
Tbxl over-expression (Supplementary Material, Fig. S2). Luci- 
ferase activity increased in response to Gata4 or Ml transfec- 
tion of C2C12 cells, confirming the ability of Gata4 and Isll 
to activate Mef2c transcription (35) also in this system 
(Fig. 5A and B). Transfection of a Tbxl expression vector 
alone did not affect luciferase activity (Fig. 5A). However, 
co-transfection of Tbxl and Gata4 showed that Tbxl inhibited 
the ability of Gata4 to transactivate the reporter (Fig. 5C). In 
contrast, Tbxl did not inhibit the Isll -induced activation of 
the reporter (Fig. 5D), suggesting that Tbxl may regulate 
Mef2c expression through a Gata4-dependent mechanism. 

It has been shown that Gata4 is up-regulated in Tbxl~'~ 
mutants (13). We analyzed Gata4 expression at E8.5 and 
E9.5 by whole-mount ISH on control and COKT;Mespl Cre/+ 
mutant embryos, and we found a marked Gata4 down- 
regulation in COET ;Mespl Cre/+ mutant hearts (Fig. 6C and 
D) compared with control embryos (Fig. 6A and B), thus con- 
firming that Tbxl can negatively affect Gata4 expression. 
Since Mef2c is regulated by Gata4 (13), one can argue that 
reduced Mef2c expression observed in our experiments is 
merely the result of Gata4 down-regulation by Tbxl. 
However, this is unlikely since we performed our luciferase 
assays on cells transfected with a Gata4 expression vector, 
which is not regulated by Tbxl. Indeed, co-transfection of 



Tbxl did not affect the amount of protein produced by the 
Gata4 expression vector (Supplementary Material, Fig. S3). 



Gata4 is required for Tbxl-induced regulation of Mef2c in 
C2C12 cells 

The previous experiments support the existence of a 
Tbxl-|Gata4^Mef2c pathway, but this is not sufficient to 
explain all our data. To test whether Tbxl regulates Mef2c 
transcription through a mechanism involving Gata4, we eval- 
uated the response of Mef2c to Tbxl in the absence of Gata4. 
C2C12 myoblasts express endogenous Gata4, so, using 
Gata4-specific siRNAs, we knocked down the mRNA and 
protein levels to ~20% of the wt level (Fig. 7 A and B). 
Gata4 knockdown was associated with down-regulation of 
Mef2c to ~80% of the wt level (Fig. 7C). Transfection of 
Tbxl (Fig. 7D) did not result in further significant reduction 
in Mef2c expression (Fig. 7C), suggesting that Gata4 is 
required for Tbxl-induced down-regulation of Mef2c. 

Tbxl occupies conserved regions of the Mef2c gene 

We performed a bioinformatic analysis to look for potential 
TBEs evolutionary conserved between the human and the 
mouse, in the Mef2c locus (see Materials and Methods). One 
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Figure 5. Tbxl inhibits Gata4-dependent activation of the MEF2 reporter construct in luciferase assays. (A and B) Gata4 and Isll are able to activate the MEF2 
reporter construct in transfected C2C12 cells (A and B) while Tbxl is not able to repress it (A). (C and D) Tbxl is able to contrasts Gata4- (C) but not Isl 1 -dependent 
activation of the MEF2 reporter construct (D). Luciferase activity was normalized to |3-galactosidase activity. Values are means + SE, n = 3,*P < 0.05. 
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Figure 6. Tbxl affects Gata4 expression in vivo. (A-D) Whole-mount ISH 
showing Gata4 down-regulation in the heart of COET; Mespl c,e/+ mutant 
embryos (C and D) compared with wt (A and B) littermates at E8.5 and 
E9.5. Arrows indicate the inflow region of the heart. 



of the putative TBEs that we have identified resides in a skel- 
etal muscle-specific control region (36) upstream of Mef2c 
(Supplementary Material, Fig. S4A). 

We also found a putative TBE in a conserved SHF-specific 
enhancer (Supplementary Material, Fig. S4B), previously 
described by Dodou et al. (35). 

Chromatin immunoprecipitation (ChIP) experiments on 
C2C12 cells revealed that endogenous Tbxl occupies both 
these sites (Fig. 8), suggesting that it may have a direct role 
in Mef2c transcriptional regulation. 

Next, we carried out luciferase assays using a reporter con- 
struct carrying the SHF enhancer (SHF-Enh-Luc; Fig. 9A) to 
test its response to Tbxl in C2C12 cells. The enhancer 
responded to transfected Isll and Gata4 but not to transfected 
Tbxl alone (Fig. 9B). In addition, as shown in the experiment 
with the desMEFl reporter, Tbxl did not inhibit Isll -induced 
activation, but did interfere with Gata4-induced activation 
(Fig. 9C and D). The mutation of the TBE site resulted in a 
significant reduction in the response to Gata4 and in a 
milder, but still significant response to Tbxl co-transfection 
(Fig. 9D). Thus, in this assay, the TBE is not completely 
required for Tbxl to carry out its negative effect upon Gata4 
activation, although this effect appears to be milder. 
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Figure 7. Gata4 is required for Mef2c repression by Tto/ in C2C12 myoblasts. (A and B) siRNA-mediated knockdown of Gata4 results in a reduction in mRNA 
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Figure 8. Endogenous Tbxl binds conserved regulatory regions in the Meflc 
locus. Small panel top-right: WB analysis of C2C12 nuclear extracts using an 
anti-Tbxl antibody demonstrating endogenous expression of the protein. 
Graph: chromatin immunoprecipitation analysis of the Meflc skeletal 
muscle and SHF TBEs performed in C2C12 myoblasts using an antibody 
against Tbxl together with rabbit IgG as the background control. qRT-PCR 
analysis of an immunoprecipitated material was performed and percent of 
input for each sample is reported in the graph. Values are means + SD, n = 3. 



DISCUSSION 

Dosage of Tbxl correlates well (though often not linearly) to 
phenotypic severity in vivo (9). With this work, we initiated 
an effort to correlate phenotypic changes with genome-wide 
transcriptional analysis, in order to identify candidate effectors 
that may be closely correlated with different phenotypic fea- 
tures. The experimental approach used here has some limita- 
tions, e.g. the whole embryo RNA analysis might smooth 
out or even hide tissue-specific changes of expression. In add- 
ition, stage-specific differences cannot be detected by our 
single-stage analysis. Nevertheless, our study was sufficient 
to identify a relatively subtle expression modulation of the 
Mef2c gene, which encodes a transcription factor critical in 
cardiac differentiation. In vivo and in vitro gain- and 
loss-of-function studies confirmed the negative correlation 
between Tbxl and Mef2c expression. This finding is consistent 
with the previously determined ability of Tbx 1 to inhibit car- 
diomyocyte differentiation in vivo and in vitro (14). Our study 
also indicates that Tbxl over-expression has inhibitory effects 
upon Mef2c expression during skeletal muscle differentiation. 

We asked how Tbxl down-regulates Mef2c expression. Our 
data suggest three, non-mutually exclusive, mechanisms 
(Fig. 10). One is direct interaction with the Mef2c gene 
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Figure 9. Response of the Mef2c SHF enhancer to Tbxl expression in luciferase assays. (A) Scheme of the Mef2c SHF enhancer- luciferase construct. The 
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obvious. Luciferase activity was normalized to (3-galactosidase activity. Values are means + SE, n = 3, *P < 0.05. 



(Fig. 10), as suggested by occupation of two TBE sites by 
Tbxl. The consequences of this interaction remain to be inves- 
tigated, but they may be chromatin changes leading to tran- 
scriptional down-regulation. Our luciferase assay with a 
mutated TBE in the SHF resulted in very limited reduction 
in the capacity of Tbxl to inhibit enhancer activity. 
However, because enhancer constructs are not properly chro- 
matinized and are out of the normal context of the endogenous 
gene, it is difficult to draw conclusions from these data. 

The second mechanism is mediated by a repressive effect of 
Tbxl upon Gata4 gene expression (Fig. 10), as seen in vivo 
(13) (data presented here). This mechanism is supported by 
tissue culture data, showing that silencing Gata4 suppresses 
the ability of Tbxl to down-regulate Mef2c in C2C12 cells. 
This effect would have negative consequences upon the ex- 
pression of MeJ2c because this gene is a target of Gata4. 
However, our desMEF2 reporter data indicate that Tbxl is 
able to inhibit reporter activation in the presence of transfected 
Gata4, and in this scenario, the transcriptional regulation of 
the Gata4 gene is irrelevant. It is more likely that Tbxl inter- 
feres negatively upon the ability of the Gata4 protein to 



regulate Mef2c (Fig. 10). The third mechanism is also compat- 
ible with the results obtained by Gata4 knockdown (Fig. 7). 
Tbx20 synergizes with Gata4 (and Isll) to activate both the 
Mef2c and Nkx2.5 SHF enhancers (37); Tbxl may compete 
with another T-box protein and inhibit this synergy. Interest- 
ingly, in our luciferase assay with the SHF enhancer, the mu- 
tation of the TBE site significantly reduced the ability of Gata4 
to activate the enhancer. This might be because of an endogen- 
ous T-box protein that synergizes with Gata4 in C2C12 cells. 
However, this is not Tbx20, as the gene encoding it is not 
expressed in these cells (our unpublished data). In addition, 
Gata4 and Tbx5 co-immunoprecipitate and cooperatively 
transactivate a luciferase atrial natriuretic factor reporter 

(38) . However, Tbxl does not co-immunoprecipitate with 
Gata4 in C2C12 cells (data not shown). Of note, 
co-transfection experiments showed the ability of Nkx2-5 
and GATA4 to transactivate the minimal Cx40 promoter, 
while Tbx5 repressed Nkx2-5/GATA4-mediated activation 

(39) , demonstrating that T-box transcription factors could ac- 
tivate or repress transcription, depending on the specific 
context. 
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Figure 10. Schematic model of Mef2c regulation by Tbxl. Our data suggest 
three, non-mutually exclusive, mechanisms: (i) direct transcriptional regula- 
tion of Mef2c by Tbxl; (ii) negative regulation of Gata4 expression and (iii) 
interference with Gata4-dependent activation of Mef2c (double-dashed line). 



To allow cells to differentiate, medium was changed to 
DMEM medium supplemented with 2% horse serum 
(Hyclone). 

We have previously described the pCDNA-Tbxl-c-myc ex- 
pression vector (11). The <ie.sMEF2-luc reporter construct (33) 
was kindly provided by Dr Eric N. Olson; the 
pCDNA3.1-Gata4-c-myc-His and pCDNA3.1-Ilsl were 
kindly provided by Dr Robert J. Schwartz. 

The SHF-Enh-Luc reporter construct was generated by 
cloning the SHF minimal enhancer of the mouse Mef2c gene 
described by Dodou et al. (35) into the pGL3- Promoter 
Vector (Promega). The mutated SHF-Enh-Luc reporter con- 
struct was generated by substitution of three nucleotides in 
the conserved TBE of the original SHF-Enh-Luc reporter con- 
struct, using the QuikChange II Site- Directed Mutagenesis Kit 
(Stratagene) according to the manufacturer's instructions. 

C2C12 cells were transiently transfected using the PolyFect 
Transfection Reagent (QIAGEN) according to the manufac- 
turer's instructions; cells were collected 24 h after transfection 
for RNA or protein extraction. 

For differentiation experiments, 24 h after transfection, 
medium was changed to differentiation medium and cells 
were collected at different days for RNA extraction. 



In summary, our results illustrate the power of multiallelic 
gene expression analysis as it identified a relatively subtle 
regulation of a gene encoding a developmentally critical tran- 
scription factor. The mechanisms by which Tbxl regulates 
Mef2c expression appear to be complex and include direct 
interaction with the gene, regulation of Gata4 expression, 
which in turn regulates Mef2c, and possibly interference 
with the Gata4 protein ability to regulate Mef2c. 

The identification of the negative regulation of MeJ2c by 
Tbxl provides an additional clue as to how Tbxl carries out 
its negative modulation of cardiomyocyte differentiation and 
provides a new link in the genetic network operating in 
cardiac progenitor cells. 

MATERIALS AND METHODS 

Mouse lines, breeding and genotyping 

All the mouse mutant alleles used in this study have been pre- 
viously reported: Tbxl~ (7), Tbxl Neo (11), Tbxl Neo2 (12), 
COET (10) and Mespl Cre (40). The COET transgene (10) 
upon Cre-mediated recombination in average expresses a 
Tbxl cDNA approximately twice as much as the level of 
one endogenous Tbxl allele (10). 

Various homozygous and compound heterozygous mutants 
were generated by heterozygous mutant mating. Embryos 
were collected at E8.5 and E9.5, considering the day of obser- 
vation of a vaginal plug to be E0.5. Mice were genotyped by 
PCR as described in the original reports. 



Cell lines, plasmids and transfections 

Mouse C2C12 myoblasts (ATCC) were cultured in Dulbecco's 
modified Eagle's medium (DMEM, Invitrogen) supplemented 
with 10% fetal bovine serum (Hyclone) at 37°C with 5% C0 2 . 



Microarray-based gene expression analysis 

Total RNA from E9.5 embryos with nine different genotypes 
(Tbxl +/+ , Tbxl Neo2/+ , Tbxl Neo/+ , Tbxl + '-, Tbxl Neo2Mo2 , 
Tbxl Neo2/Neo , Tbxl Neo2/ -, Tbxl Neo/Neo and Tbxr'~) was iso- 
lated using TRIzol (Invitrogen) and cleaned using the 
RNeasy cleanup Kit (QIAGEN) according to the manufac- 
turer's instructions. Procedures for cRNA preparation and 
GeneChip processing were performed at Boston University 
Microarray Resource Facility exactly as described in the Affy- 
metrix GeneChip Expression Analysis Technical Manual 
(Affymetrix, Santa Clara, CA, USA, current version available 
at www.affymetrix.com). Briefly, total RNA integrity was 
verified using RNA 6000 Nano Assay RNA chips run in 
Agilent 2100 Bioanalyzer (Agilent Technologies). Total 
RNA was reverse transcribed using One Cycle cDNA Synthe- 
sis kit (Affymetrix) and the obtained double-stranded cDNA 
was purified with GeneChip Sample Cleanup Module (Affy- 
metrix). cDNA was used as a template for in vitro transcrip- 
tion using GeneChip in vitro transcription (IVT) Labeling 
Kit (Affymetrix). The obtained biotin-labeled cRNA was puri- 
fied using GeneChip Sample Cleanup Module (Affymetrix), 
fragmented and hybridized (10 |xg) to the GeneChip Mouse 
430 2.0 microarrays (Affymetrix) for 16 h in GeneChip Hy- 
bridization oven 640 at 45°C. For each genotype, we used 
two biological replicates, each hybridized to one array, with 
a total of 1 8 arrays. IVT and cRNA fragmentation quality con- 
trols were carried out by running an mRNA Nano assay in the 
Agilent 2100 Bioanalyzer. The hybridized samples were 
washed and stained using Affymetrix fluidics station 450 
with streptavidin-R-phycoerythrin (SAPE) and the signal 
was amplified using a biotinylated goat anti-streptavidin anti- 
body followed by another SAPE staining (Hybridization, 
Washing and Staining Kit, Affymetrix). Microarrays were im- 
mediately scanned using Affymetrix GeneArray Scanner 3000 
7G Plus (Affymetrix). The data were analyzed with the 
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Bioconductor package 'Affy' using Robust Multiarray 
Average method for normalizing and summarizing probe 
level intensity measurements. 



Statistical analysis of microarray data 

We developed a strategy for statistical comparison of chip data 
aimed at identifying genes that responded to variations of Tbxl 
dosage. To discover which genes responded significantly, we 
computed the log ratio of all dosage points with respect to 
the control (wt) point for all genes. After this, we used a 
two-step filtering: in the first step, we made the assumption 
that if a gene is not responding to the dosage variation, then 
the area under the gene profile curve should be close to 
zero. We approximated this area by summing the square of 
the relative expression levels of all genes. Under the null hy- 
pothesis that a gene is not responding in the dosage series, this 
area should be close to zero, otherwise it should be higher. We 
approximated the area distribution with an exponential prob- 
ability density function with the mean equal to the average 
of the area of all transcripts in the microarray. By using a 
P-value of 0.05, we selected a list of 2230 transcripts, includ- 
ing Tbxl (Supplementary Material, Table SI). 

In the second step, we performed a more stringent test based 
on a Bayesian analysis of signals (41). This fully Bayesian ap- 
proach allowed us to automatically identify and rank 497 dif- 
ferentially significant genes, whose expression is sensitive to 
Tbxl dosage (Supplementary Material, Table S2). 



GO and cluster analyses 

GO analysis (http://david.abcc.ncifcrf.gov/) was performed to 
identify specific biological pathways enriched in the list of 498 
differentially expressed genes. The most significant terms 
included regulation of metabolic processes and regulation of 
gene expression (false discovery rate < 0.05). 

Hierarchical clustering was performed on the list of 498 dif- 
ferentially expressed genes using the MATLAB software (http 
://www. mathworks.it/help/toolbox/stats/bq_679x-3. htm 
l#bq_6_ia). Four main clusters were identified and, for each 
cluster, we calculated the average profile and the associated 
standard deviation (SD). 



TBE search 

Upstream and intronic genomic sequences for human and 
mouse Mef2c were obtained from both the Ensembl and 
UCSC Genome Browser web sites and were aligned using 
the multisequence local alignment tool MULAN (http://mula 
n. dcode.org/). No specific Tbxl -binding consensus sequence 
has been identified so far; however, Tbxl can activate, in 
vitro, the FgflO promoter, through a conserved Tbx5-binding 
site (11). Multi-TF (http://multitf.dcode.org) transcription 
factor-binding site analysis was performed to find evolutionary 
conserved Tbx5-binding elements, in the aligned sequences. 
With this approach, we uncovered several TBEs in the 
Mef2c locus. 



Quantitative real-time PCR 

Total RNA was extracted from whole E9.5 wt and mutant 
embryos using TRIzol (Invitrogen) and from C2C12 using 
the RNeasy Mini Kit (QIAGEN). One microgram of total 
RNA was reverse transcribed using Superscript III (Invitro- 
gen) according to the manufacturer's instructions, in a total 
volume of 20 the resulting cDNA was diluted 1:10. 
qRT-PCR was performed using SYBR Green for the detec- 
tion of fluorescence during amplification. Each amplification 
reaction contained 25 |xl of Platinum SYBR Green qPCR 
SuperMix-UDG (Invitrogen), 0.2 |xg each of forward and 
reverse primers and 1 jjlI of diluted cDNA. All primers were 
designed with annealing temperatures of 58-62°C. PCR con- 
ditions were: 50°C for 2 min, 95°C for lOmin, 40 cycles of 
95 °C for 15 s, 60°C for 1 min and followed by a dissociation 
stage of 95°C for 15 s and 60°C for 15 s. All the genes were 
amplified on a 7900 HT Fast Real-Time PCR System 
(Applied Biosy stems). PCR amplifications were performed 
in triplicate and at least three separate analyses were carried 
out for each gene. Melting curve analysis was performed 
after each run to check for the presence of non-specific PCR 
products and primer dimers. The qRT-PCR results, recorded 
as threshold cycle numbers (C t ) were normalized using 
GAPDH as an internal control and analyzed using the com- 
parative C t method (2~ AACl method) (42,43). 



In situ hybridization 

Digoxigenin-labeled RNA probes were prepared by standard 
methods (Roche). Mef2c and Gata4 probes were kindly pro- 
vided by Dr Brian L. Black. Embryos were dissected in 
di-ethylpolycarbonate (Sigma Aldrich)-treated phosphate- 
buffered saline (PBS) and fixed overnight in 4% paraformalde- 
hyde-PBS (Sigma Aldrich) at 4°C. 

For ISH on cryosections, embryos were cryoprotected with 
increasing concentrations of sucrose (Sigma Aldrich), treated 
with a mixture of 1 part 30% sucrose and 1 part OCT for at 
least 2 h at 4°C, embedded in OCT and cut to generate 10- 
12 |jLm sections. After post-fixation in 4% paraformalde- 
hyde-PBS, sections were incubated in triethanolamine 
buffer (0.2% acetic acid and 1% triethanolamine, Sigma 
Aldrich) containing acetic anhydride and washed twice in 
PBS for 5 min. The sections were then prehybridized for 1 h 
at 70°C in the hybridization mixture (50% formamide from 
Sigma Aldrich, 5 x saline-sodium citrate (SSC), salmon 
sperm DNA 40 jig/ml and 25 mg/ml yeast tRNA from Invitro- 
gen). The probes were denatured for 5 min at 80°C and added 
to the hybridization mix (400 ng/ml). The hybridization reac- 
tion was carried out over night at 70°C; prehybridization and 
hybridization were performed in a box saturated with a 5 x 
SSC-50% formamide solution to avoid evaporation. After in- 
cubation, the sections were washed twice with MABT (50 min 
maleic acid, pH 7.5, 250 mM NaCl and 0.1% Tween-20, Sigma 
Aldrich), treated with MABT containing 10% sheep serum for 
1 h at room temperature and then incubated with alkaline 
phosphatase-coupled anti-digoxigenin antibody (Roche) 
diluted 1:2000 in MABT containing 10% sheep serum, over- 
night at 4°C. The day after, sections were washed with 
MABT and equilibrated for 5 min in Buffer 3 (Tris-HCl 
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100 mM, NaCl 100 mM and MgCl 2 50 him, pH 9.5); color de- 
velopment was performed at room temperature (time depend- 
ing on the amount of transcripts to be detected) in Buffer 3 
containing nitro blue tetrazolium and 5-bromo-4-chloro-3- 
indolyl phosphate (Roche). Staining was stopped by a wash 
in PBS containing 0.1% Tween-20; sections were rehydrated 
for 5 min in deionized water and then dehydrated through suc- 
cessive baths of EtOH (70, 95 and 100%) and xylol, and 
mounted in Eukitt resin (O. Kindler GmbH & Co.). 

Whole-mount ISH was performed according to the previ- 
ously published methods (44). At least three somite staged 
wt and mutant embryos were analyzed with each probe. 

Western blotting 

Nuclear proteins were extracted as described previously (18). 
The following primary antibodies were used: rabbit anti-Tbxl 
(Abeam, 1:1000), goat anti-Mef2c (Santa Cruz Biotechnology, 
1:1000), rabbit anti-Gata4 (Santa Cruz Biotechnology, 
1:1000), monoclonal anti-c-myc (Sigma- Aldrich 1:1000) and 
monoclonal anti-fJ-Actin (Sigma- Aldrich 1:10 000). The fol- 
lowing secondary antibodies were used: horseradish 
peroxidase-conjugated (HRP) anti-rabbit and anti-mouse (GE 
Healthcare, 1:10 000) and HRP-conjugated anti-goat (Santa 
Cruz Biotechnology, 1:5000). The HRP-derived signal was 
detected using the Amersham ECL and ECL Plus Western 
Blotting Detection Reagents (GE Healthcare). 

Luciferase assay 

For luciferase assays, mouse C2C12 myoblasts were plated in 
24-wells and transfected with PolyFect Transfection Reagent 
(QIAGEN) according to the manufacturer's instructions. The 
afesMEF2-luc reporter construct and, for single experiments, 
the pCDNA-Tbxl-c-myc, pCDNA3.1-Gata4-c-myc-His and 
pCDNA3.1-Ilsl plasmids were always co-transfected with a 
pCMV-(3-Gal expression vector (Clontech); at 24 h after trans- 
fection, cell extracts were prepared and activities of 
(3-galactosidase and Firefly luciferase were measured using 
the Luciferase Assay System (Promega) and the Beta-Glo 
Assay System (Promega), respectively, on a GLOMAX 96 
microplate luminometer (Promega). Relative luciferase activ- 
ity (Firefly luciferase for reporter and (3-galactosidase for nor- 
malization of transfection efficiency) was measured following 
manufacturer's instructions (Promega). The data represent the 
means and SDs of at least three independent transfections. 

RNA interference 

Gata4 interference in C2C12 myoblasts was achieved by using 
a set of four different siRNAs specifically directed against the 
mouse Gata4 locus (ON-TARGET plus SMART Pool, 
Thermo Scientific). As negative control, we used a pool of 
siRNAs that virtually targeted no mouse genes (ON-TARGET 
plus Negative Control, Thermo Scientific). Cells were plated 
in six wells and transfected with 50 nM/well of siRNAs 
using Lipofectamine 2000 (Invitrogen) according to the man- 
ufacturer's instructions. The Gata4 mRNA level was mea- 
sured by qRT-PCR and the Gata4 protein level by WB as 
described above. 



ChIP assay 

ChIP assay was performed as described previously (45). Son- 
ication was performed for 10 cycles of 15 s on and 60 s off 
using a Soniprep 150 (Measuring and Scientific Equipment). 
Immunoprecipitation was performed overnight at 4°C with 
an antibody against Tbxl (Abeam, 4/100 |xg of chromatin); 
normal rabbit immunoglobulin G (IgG) was used as the back- 
ground control. 

Results were quantified by qRT-PCR using the primers 
listed in Supplementary Material, Table S3, and the enrich- 
ment of DNA was calculated in terms of percent of input, as 
described previously (46). PCR amplifications were performed 
in triplicate from three independent experiments. 

Statistical analysis 

All data were expressed as means + standard error or SD from 
three independent experiments. Differences between groups 
were examined for statistical analysis using a two-tailed Stu- 
dent's f-test. Values of P < 0.05 were considered significant. 

SUPPLEMENTARY MATERIAL 

Supplementary Material is available at HMG online. 
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